DECOHERENCE-INSENSITIVE QUANTUM 
COMMUNICATION BY OPTIMAL C*-ENCODING 



BERNHARD G. BODMANN 1 , DAVID W. KRIBS 2 ' 3 , AND VERN I. PAULSEN 4 



Abstract. The central issue in this article is to transmit a quantum 
state in such a way that after some decoherence occurs, most of the in- 
formation can be restored by a suitable decoding operation. For this 
purpose, we incorporate redundancy by mapping a given initial quantum 
state to a messenger state on a larger-dimensional Hilbert space via a 
C*-algebra embedding. Our noise model for the transmission is a phase 
damping channel which admits a noiseless or decoherence-free subspace 
or subsystem. More precisely, the transmission channel is obtained from 
convex combinations of a set of lowest rank yes/no measurements that 
leave a component of the messenger state unchanged. The objective of 
our encoding is to distribute quantum information optimally across the 
noise-susceptible component of the transmission when the noiseless com- 
ponent is not large enough to contain all the quantum information to be 
transmitted. We derive simple geometric conditions for optimal encoding 
and construct examples. 



1. Introduction 

The use of quantum systems as a medium to carry information enables 
algorithms and communication techniques not possible in the purely classical 
setting. While there is considerable excitement about these new capabilities, 
their implementation poses severe theoretical and experimental challenges. 
One of the major impediments to large-scale experimental advances is the 
decoherence of quantum systems as they evolve in time. Indeed, finding 
efficient methods for "quantum error correction" is of central importance in 
the fields of quantum computation, cryptography and communication. 

A fundamental, "passive" technique for quantum error correction grew out 
of the recognition that certain error models contain symmetries that can be 
used to hide encoded qubits from the overall noise of the model [1-7]. This 
scheme relies on the identification of what are now known as noiseless or 
decoherence-free subspaces or subsystems. (For convenience we will just use 
the noiseless terminology.) A number of experimental advances have been 
based on this approach [8-11]. Noiseless subsystems and subspaces have 
recently received wider attention, such as in quantum information processing 
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and cryptography [12-14], and in quantum cosmology and gravity [15-17]. 
Moreover, a structure theory is now emerging [18-23] which shows, among 
other things, how such subsystems and subspaces can be identified [22,23]. 

The basic problem we address in this paper is, given a certain class of 
decoherence error models, how to distribute quantum information optimally 
across the noise-susceptible component of the transmission when the noiseless 
component is not large enough to encode all the information to be transmit- 
ted. Our analysis is motivated by recent work [24-28] on the construction 
and use of frames to incorporate redundancy in classical signal transmission. 

For the purposes of this paper, a quantum communication framework is 
given by triples from a collection of quantum channels. Each triple contains 
the maps for encoding, transmission and reconstruction of quantum states. 

We exclusively consider encoding maps that are trace-preserving C*-alge- 
bra embeddings in order to preserve the structure of the set of encoded quan- 
tum states. We call such maps C*-encodings. Consequently, the encoding 
and the subsequent reconstruction are quantum channels, that is, completely 
positive and trace-preserving (or at least non-increasing for the trace). The 
encoding channel incorporates redundancy by mapping a quantum state on 
a given Hilbert space to a "messenger" state on a larger-dimensional Hilbert 
space. 

It is common in the context of quantum error correction to model noise in 
transmissions by a so-called interaction algebra A and search for states that 
are in the commutant of A. The interaction algebra we consider is given by 
the commutative set of operators A = {D ® I : D diagonal} on a Hilbert 
space C m <S> C l , admitting a noiseless subsystem in the second component. 
However, we do not consider the entire commutative interaction algebra A, 
but rather focus on a "minimal" set of error models for the transmission. 
This set contains all nontrivial lowest-rank projective measurements and 
their convex combinations. Such errors model generalized phase damping, 
that is, decoherence due to minimal interaction with a macroscopic system 
that performs with certain probabilities one among several possible yes/no 
measurements on our system. The commutativity means that the outcomes 
of the measurements obey the rules of classical Boolean logic. Apart from 
naturally occurring decoherence due to interactions with an environment, 
the motivation for considering this type of error is that such measurements 
may be used to remove corrupted information in a quantum state after the 
occurrence of a quantum error has been detected. 

Minimal projections in the interaction algebra A are tensor products 
of diagonal matrix units {Ejj}™ =1 and the identity, henceforth denoted as 
Qj = Ejj ® /. Accordingly, a transmission channel which performs a mini- 
mal projective measurement of a messenger state M is : M i— > QjMQj + 
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Qj-MQj-, where Q\ complements the orthogonal projection Q~ to the iden- 

J J J ' ' "* 

tity, Qj + Qf = I. We denote the phase damping channels obtained by 
convex combinations of such measurements as 

m m 

(1) Q (1) = {Y^pMM) I all Pj > 0, J> = l} • 

In the following, we shall use B(C d ) to denote the algebra to be en- 
coded. This algebra is embedded via the encoding map $ into the higher- 
dimensional algebra B(C m ®C l ) = B(C m )®B(C l ). The subsequent decoding 
is a "blind" reversal of the encoding operation. 

For our purposes, an optimal C*-encoding <£> minimizes sup^- e(<3>, £), the 
supremum of the Hilbert-Schmidt norms for the reconstruction error occur- 
ring among all quantum states and all noisy transmission channels 8 in the 
convex set While other measures for distortion are also used in quan- 

tum information, we have chosen the Hilbert-Schmidt norm because it is 
most closely related to the geometry of trace-preserving C*-algebra embed- 
dings. 

Our main result is that if I < d, a C*-encoding <3> satisfies 

sup e(*S)/d>l 2 - 2d l ml > , d/ml e [0,1/2]; 

S) I (8mld-2m 2 l 2 -Ad 2 )~ 1 /\ d/ml € (1/2, 1] . 

In addition, we characterize the case of equality as that for which the encod- 
ing map is constructed from what we call a family of "uniformly weighted 
projections resolving the identity". In many cases we show how such fam- 
ilies can be constructed as an extension of recent work on classical signal 
encoding problems [25]. However, we leave the problem of determining the 
possible existence of these families in complete generality for investigation 
elsewhere. Finally, for all C*-encodings that are optimal for QP~\ we derive 
an inequality for transmission channels with next-to-minimal decoherence 
and characterize the case of equality as that for which all pair sums of the 
underlying projections have equal operator norm. This geometric condition 
on the ranges of the projections generalizes the notion of equiangular lines. 

This paper is organized as follows: In Section 2 we describe the concept 
of a quantum communication framework and examine the structure of this 
framework in the presence of a noiseless subsystem. In Section 3, we obtain 
precise estimates on the reconstruction error due to a noisy transmission and 
use this to characterize optimal encoding in terms of uniformly weighted pro- 
jective resolutions of the identity. Section 4 treats optimality for secondary 
decoherence errors in the uniform encoding case. Finally, we comment on a 
generalization that includes the noiseless subspace case, with similar results 
on the reconstruction errors and optimal encoding schemes. 
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2. Quantum Communication Framework 

We recall that for the purposes of this paper, a quantum communication 
framework is given by three quantum channels: encoding, transmission, and 
decoding. Hereby, a quantum channel is understood to be a completely pos- 
itive, trace preserving (or at least trace non-increasing) linear map between 
C*-algebras. We point the reader to [29-31] as entrance points into the lit- 
erature on completely positive maps and their use in quantum computation 
and information. 

2.1. C*-Encoding and Decoding. Apart from being quantum channels, 
we require that these maps preserve the algebraic structure of the quantum 
states to be encoded. 

Definition 2.1. Let TL and JC be finite- dimensional Hilbert spaces. We say 
that a map $ : B{7i) — > B(JC) is a C* -encoding of BiTi) if $ is a trace- 
preserving *-monomorphism. The associated decoding map is the Hilbert- 
Schmidt dual $* : B(K) -> B(H). 

Remarks 2.2. Algebra embeddings of the type introduced here have also been 
used for a unified treatment of a hybrid quantum/classical memory [32]. 

In the context of our quantum communication framework, the most impor- 
tant property of the encoding is that it preserves the Hilbert-Schmidt inner 
product. Thus $* is a left inverse of $, and encoding followed by decoding 
is the identity map, = id. Moreover, the dual $* is completely positive 
and non-increasing for the trace. 

Another consequence of $ being a Hilbert-Schmidt isometry is that, in- 
stead of performing a measurement after decoding a transmitted state, one 
may simply encode the measurement and apply it directly to the transmitted 
state. 

Furthermore, the benefit of encoding with a C*-algebra embedding is that 
a quantum algorithm performed on the initial state, implemented by conju- 
gation with a unitary operator, could be performed either prior to encoding 
or afterwards. 

2.2. Error Model for Transmission. The transmission of an encoded 
state M G B(JC) is described by a quantum channel. When dim/C < oo, such 
a map has a representation of the form E(M) = ^2 a E a ME* for a finite set 
of operators {E a } C B(7i). These operators {E a } are said to form an "error 
model" for the channel. It is the impact of such error operators that must 
be overcome to safely transmit quantum information through the channel 8. 

As discussed in the introduction (see Eq. (1)), here we consider phase 
damping channels given by convex combinations of commuting projective 
measurements. 
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Definition 2.3. We call a quantum channel £ a generalized phase- damping 
channel on B{K) if there is a commutative set of projections {Qj}Y=i C B(JC) 
and a probability vector p G M™, J2jPj = 1, such that for M G B(JC), 

m 

£{M) = J^PAQjMQi + QfMQf) , 
j'=i 

where each Qj- is the projection on the orthogonal complement of the range 

ofQj. 

We define the overall reconstruction error of the framework as follows. 

Definition 2.4. Let $ : B(H) — > -B(/C) fre a C* -encoding as in Defini- 
tion 2.1, and let £ : B{K) — > £>(/C) fre a completely positive, trace preserving 
map on the Hilbert space K. We then define the reconstruction error caused 
by £ in the course of transmitting a state W G B(H) as 

Y = W - ($* o£o$)(W) 

and the maximum of the Hilbert- Schmidt norm ofY over the set of all states 
{W > 0, trW = 1} as the worst-case error norm 

e($,£) = max(tr[FF*]) 1/2 . 
w 

Remark 2.5. The Hilbert-Schmidt norm is not directly motivated by informa- 
tion-theoretic considerations. However, in the situation considered here, it 
is a natural way to quantify the impact of an error £ in the transmission 
with $* o £ o $. In fact, the Hilbert-Schmidt norm is non-increasing under 
any generalized phase-damping channel £. This follows from the convexity 
of the norm and from the case of a single projection Q and a state M, 

tr[(QMQ + Q ± MQ ± ) 2 ] = ti[QMQ 2 MQ + Q ± M(Q ± ) 2 MQ ± ] 

< tr[QM 2 Q + Q ± M 2 Q ± ] = tr[M 2 } . 

Since this norm is also non-increasing when $* is applied to £(M), we can 
use it to measure the amount of corruption caused by the error £ and the 
subsequent decoding. When calculating the Hilbert-Schmidt norm of the 
maximal reconstruction error, we thus bound the maximal extent of corrup- 
tion for any state. 

2.3. Encoding in the Presence of Noiseless Subsystems. Consider a 
Hilbert space /C for a quantum system and a channel £ : B{K,) — > B(JC). 
Suppose we have a decomposition K = (/C (1) <g>/C (2) )©V, so V = (/C^^/C^)- 1 , 
and for all states W w and W {2) on /C (1) and /C (2) there is a state S {1) on /C (1) 
with 

£(W {1) ® W {2) © 0) = S {1) © W {2) © . 
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Then the second component of the tensor product is said to be a noiseless 
subsystem (or noiseless subspace in the case dim/C* 1 ) = 1) for £. The ter- 
minology "decoherence-free" subspace and subsystem is also quite common. 
This definition ensures that quantum information encoded in the second 
component of the tensor product is unaffected by the noise of the channel £. 

For succinctness, we only treat the case V = {0} here. Remark 4.5 sketches 
how the results generalize to the case dimV > 0. 

For our transmission channels, the second component of the tensor product 
decomposition C N = C m <8> used for the transmission is noiseless, whereas 
the first component may lose its off-diagonal entries. If / > d, then we could 
simply encode W as M = —I (&W, which is unaffected by any decoherence 
in the first component. So let us assume I < d, which means the dimension 
of the second component is too small to contain all the information given by 
an arbitrary state W on C d . 

2.3.1. C* -encodings and Positive Operator Valued Measures. Let us describe 
a useful concrete description of C*-encodings in terms of associated factor- 
izations of POVMs. 

Lemma 2.6. Given a C* -encoding $ : B(C d ) -> B(C m ®C l ), then there 
exists a set of so-called coordinate operators {Vj : C d — > C'} resolving the 
identity Y^j=i = I suc ^ that $ has the form 

$ : B[C d ) -> B(C m <g> C l ) = B(@f =l C l ) 

F » G = (G t3 )™ =1 ,G tJ = VWV;. 

Moreover, the action of the decoding map is given by 

m 

(2) G^$*(G) = J2 V * G v V i- 

Proof. As $ is a trace-preserving C*-monomorphism of a simple, finite- 
dimensional C*-algebra B(C d ), it follows from the representation theory for 
such algebras [33] that there is an isometry V : C d -> C rn ® = ®f =1 C l , 
V*V = J, which for all F e B(C d ) gives 

(3) = VFV* . 

We may view V as a column matrix of coordinate operators Vj : C d — > C', 
1 < j < m, 

(4) V = {V X V 2 ■■■ V m y : C d — . ef =1 C l ; / ^ (VJ V 2 f ■ ■ ■ V m f) 1 . 

The linearity of V ensures that each Vj is linear. Also note that rank(Vj) < /, 
1 < j < m. 
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Thus, the action of $ is given by G = $(F) = VFV*. We may view G as 
an operator on C m <E> C expressed in block matrix form 

(5) C7 = (cy™ =1 with g„ = v^w/ . 

Now the Hilbert-Schmidt dual of $ can be identified as 

m 
i,3=l 

□ 

Remarks 2.7. The operators {V^-} define a "positive operator valued mea- 
sure" (POVM), {A,-}, where A,- = V*Vj and 

m 

(6) I = V*V = Y,Ay 

3=1 

Conversely, if we start with a POVM on C d given by {A,} such that Aj > 0, 
each Aj has rank at most I, and Y^j=i A? = I ■> u P on factoring Aj = V*Vj 
one can define V as in Eq. (4) and a C*-encoding as in Eq. (3). 

Hence the C*-encodings of Definition 2.1 are completely characterized by 
factorizations of POVMs in this way, and we shall use this concrete form in 
the analysis below. A special case is of central importance in this paper. 

Definition 2.8. Let H be a d- dimensional complex Hilbert space and let 
{Aj}Y=\ C B{7i) be a POVM onTi. If each operator is given by Aj = kjPj 
with a non-negative coefficient kj and a projection operator Pj = PjPj, then 
we call {Aj} a set of weighted projections resolving the identity. If kj = k = 
dj Y^j=i t r Pj f or a M j {1)2,..., m}, we call {kPj}™ =l a set of uniformly 
weighted projections resolving the identity. 

Remarks 2.9. (i) Suppose all operators in a POVM {Aj}™ =l have a maximal 
rank /, and we denote Pj the projection onto the range of each Aj, and 
abbreviate the operator norms kj = \\Aj\\. Then from kjPj > Aj and from 
taking traces on both sides of Eq. (6), we have 

m m 

(7) lJ2kj>J2 tTA i = d - 

3=1 3=1 

Equality holds in this inequality if and only if Aj = kjPj, and Pj has rank / 
for all 1 < j < m. 

(ii) Weighted projections giving a resolution of the identity and their ro- 
bustness against perturbations have been investigated by Casazza and Ku- 
tyniok in the context of frames and distributed processing, where they are 
known as Parseval frames for subspaces [34] or fusion frames [35]. The C*- 
encoding derived from this special case of POVMs is of central importance 
in this paper. 



8 



B. G. BODMANN, D. W. KRIBS, V. I. PAULSEN 



(iii) If X^=it r Pj > d, then the ranges of the P/s are not mutually or- 
thogonal. However, this does not imply that the ranges of the P/s have to 
intersect nontrivially, see the example below. 

(iv) Choosing orthonormal bases for the range of each Pj and scaling the 
basis vectors spanning the range of each Pj by a factor ^fkj would yield a 
Parseval frame of N — J2™ =1 tr Pj vectors for the space C d . If the weights 
are uniform, so is the frame. 

(v) Given a set of weighted projections resolving the identity such that 
Sj=i t r Pj > d, we can replace each Pj by Pj — I — Pj and each coefficient 
kj by k'j = kjj (X^=i — 1) in order to arrive at a complementary resolution 
of the identity, Y^j=i k'jPj = d- 

2.3.2. Examples. We sketch the construction of uniformly weighted projec- 
tions resolving the identity. 

Example 2.10. It is fairly straightforward to construct examples of weighted 
projections resolving the identity for which all {kjPj}'j! =1 are rank-one and 
m> d. To this end, we take a rank-rf orthogonal projection G = G*G on C m 
and factor it as G = VV*, where V is an isometry V : C d — > C m . Now denote 
the m column vectors of the matrix representation for V* in the standard 
basis as {fj}^Li, choose each Pj to be the projection on the one-dimensional 
subspace of C d spanned by fj, and let the set of weights be {kj = ||/j|| 2 }. It 
is then straightforward to verify that Y^j=i kjPj = d- 

Holmes and Paulsen [25, Remark 1.1] show how to make the weights uni- 
form by rotating amongst the frame vectors. 

The next example describes how rank-one projections may be used to 
construct higher-rank weighted projections which resolve the identity. 

Example 2.11. If the Hilbert space is a tensor product C d = C l (E)C q , so d = Iq, 
and there is a set of weighted rank-one projections {kjYlj} 1 JL 1 resolving the 
identity on the second component C g , then {kjPj = kjl ® H,-} is a set of 
weighted rank-/ projections resolving the identity on C d . 

If all weights are uniform in the resolution of the identity for the second 
component, the same holds for the tensor product. 

However, there are examples of weighted projections realized on tensor 
product spaces that are not of the type {/ <S> H,}. 

Example 2.12. We now describe a set {kP\, kP m } of m uniformly weighted 
projections resolving the identity in the case that d = 21 is even. 

This example has the property that rk(Pj) = I for every j and if i ^ j, 
then the ranges of Pi and Pj intersect only in the zero vector. 
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To establish the claimed properties, we identify in the usual way B(C d ) = 
B(C l © C') = P(C 2 ) <g> B(C l ); that is, as 2 x 2 block matrices, each of whose 
entries is an / x I matrix. Let 

(8) P=(JJ) and ^=(g-| 

for all 1 < j < m, where each Cj and Sj is an / x / diagonal matrix of the 
form 

i i 

(9) Cj = ^cos(^j)^, Sj = ^2sm( y 9 iJ )E ii , 

i=i i=i 

where 9 i:j = n ^~^ , 1 < j < m,l < i < I, and {-Eji}- =1 are the diagonal 
matrix units. Using these matrices, we specify our set of m projections of 
rank /, by setting 

do) Pj=VJFUi= (j? 3i -Cj*,y 

for all 1 < j < m, so that rk(Pj) = I for all j G {1,2,... m}. 

Note that in this case, k = d/ J2j^ T (Pj) — d/(ml) = 2/1. To verify that 
T,?=i Pi = f h, observe that £; =1 cos 2 (^) = f , while ^ m =1 cos(^) 
sin( 7r( -^^ ) = 0. Finally, to verify that the ranges of a pair Pj and P,, i 7^ j, 
do not intersect nontrivially, we show ||Pj + P,-|| < 2. Since the operator 
norm is invariant under conjugation with unitaries, we observe 

(11) || Pi + Pj\\ = ||P + UiU*PUjU*\\. 

(C —S\ 
g (j J , where 

(12) C7 = cos( 7r( '~ j) )/ and S = sin( ^ ~ j) )/. 

m m 

Thus, it will be enough to verify that the matrix 

(13) p + u^PUjU* = ( J _ + Jjf 

has norm strictly less than 2. However, this last matrix decomposes into 
a direct sum of 2 x 2 Hermitian matrices of the same form and one can 
check that the eigenvalues of such a matrix are, 1 ± cos(0), with 9 = ■ 
Thus, for i ^ j, 1 < i, j < m, the angle 9 is not an integer multiple of ir, so 
the eigenvalues of these Hermitian matrices, and hence their norms, will be 
bounded above by 2. Indeed, a tight upper bound on the norm is l + cos(— ). 
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3. Uniformly weighted projections as optimal encoders for 

minimal decoherence 

The main theorem of this section characterizes optimal encoding in our 
communication framework. We prepare the result with three lemmas. 

Lemma 3.1. Given a nonnegative operator A with eigenvalues {a r } C [0, 1] 
and at least one nonzero a r , then a vector p with p r > 0, J2 r p r = 1 maxi- 
mizes 

J^((l - a r )a s + (1 - a s )a r ) 2 p r p s 

r,s 

if and only if at most two entries of p are nonzero. 

Proof. This is a standard Lagrange-multiplier argument. Suppose p is a 
maximizing vector and there is an index r such that p r — 1, then the assertion 
is true. On the other hand, assume all p r < 1. Then there exists a A e R 
such that for any r with p r ^ 0, 

- a r)oi s + (1 - a s )a r ) 2 p s = A . 

s 

Since the left hand side is quadratic in a r , there can be at most two different 
values for a r that satisfy this equation. Consequently, p r can be nonzero for 
at most two r. □ 

Lemma 3.2. Let A be a nonnegative operator with eigenvalues {a r } C [0, 1], 
the maximum of which is given by k = \\A\\. Let the function g : [0, k] 2 — > M. 
be defined by g(x,y) = ((1 — kx)y + (1 — ky)x) 2 . Then the maximum of 

G(p, x, y) = p 2 g(x, x) + 2p(l - p)g(x, y) + (1 - p) 2 g(y, y) 

over p G [0, 1] and {(x,y) : k > x > y > 0} is achieved for some p when 
x = k and y — 0. 

Proof We note that the auxiliary function 

h(x, y) = x + y — 2xy 

is harmonic in D = [0, k} 2 and piecewise linear when restricted to the bound- 
ary <90. Denote 

H(p, x, y) = p 2 h(x, x) + 2p(l - p)h(x, y) + (1 - p) 2 h(y, y) . 

Let r] Xi y be the so-called harmonic measure on <9B for the point (x,y). That 
is, for any harmonic function u that is continuous on D, it gives u(x,y) = 
J do udr] Xt y. Using a suitable convex combination of harmonic measures for 
the points (x, x), (x, y), (y, x) and (y, y), we construct p p , x , y supported on the 
boundary <9B such that f m h dp PjXjV = H(p, x, y). Since h is piecewise linear 
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on dU), this measure can be replaced by weights Wn, w i2 , w 2 ±, w 22 > at the 
corners. These weights give 

w n h(K, k) + w 12 h(K, 0) + w 2 ih(0, k) + w 22 h(0, 0) = H(p, x, y) 

and can be determined as 

wn = -^(px + (1 -p)yf , 
w 22 = \ (p(k - x) + (1 - p)(k - y)) 2 , 



and 



w 12 = w 2 i = \\p 2 x(n - x) +p(l -p){y{n - x) + x(k - y)) 



+ (1 - pfy{n - y)] . 
It can be verified that wi 2 = ^/wuw 22 and that wu + 2wi 2 + w 22 = 1. 

We observe that g = h 2 is sub-harmonic. Consequently, Jensen's inequal- 
ity yields 

G(p, x, y) < / g dfi PjXty 

J dO 

and since the restriction of g to each side of <9D is convex, applying Jensen's 
inequality again gives us 

G(p, x, y) < w u g(K, k) + w 12 g(n, 0) + w 21 g(0, k) + w 22 g(0, 0) 
= G(y/w^[, k,0) . 

Thus, G assumes its maximum when the point (x, y) is on the boundary of 
the square D. □ 

Lemma 3.3. Let $ : B{C d ) -> B(C m ) <g> B(C l ) be a C*-encoding and let 
{Vj}" l =i be the associated coordinate operators of rank I < d as defined in 
Lemma 2.6. Let Q n = E nn <g> /, for 1 < n < m, be the orthogonal projection 
onto the nth copy of C' in C m <S> = ©"i 1 C' ; and denote £ n the projective 
measurement applied to any state M G ^(©^C'), 

S n (M) = Q n MQ n + QiMQi. 

Then the maximum of the Hilbert- Schmidt norm (trfyy*]) 1 / 2 of the recon- 
struction error Y = W — $*(£ n ($(W))) over all states W only depends on 
the operator norm k n = ||V^V^||. Specifically, we have 

,1/2 

and 

(15) e($,g n ) = -== kn when k n > 1/2. 

y/-2 + 8k n - Ak 2 n 



(14) e($,£ n ) = max(tr[yy*]) 1/2 = 2k n (l - k n ) when k n < 1/2, 
w 
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Proof. We fix n and abbreviate A — A n — V*V n . After decoding, the result- 
ing error in the reconstructed state is 

Y = (I - A)WA + AW (I - A) . 

The square of the Hilbert-Schmidt norm is 

tr[yy*] = tr[((J - A)WA + AW {I - A)) 2 } 

= 2tr[((7 - A) AW) 2 + (/ - A) 2 WA 2 W] . 

In the eigenbasis of the operator A, we can express this as 

tr[yy*] =2£((1" a r )a r (l - a s )a s \W r , s \ 2 + (1 - a r ) 2 a 2 \W r , s \ 2 ) , 

r,s 

where the matrix with entries (W r>s ) represents the state W in the eigenbasis 
of A and {a r } are the corresponding eigenvalues of A. 

We notice that if we assume a r < 1 for all r, the expression for tr[YY*] 
is just the square of a weighted £ 2 -norm on the set of all states, and thus 
convex in W. Moreover, the set of all states is convex, and therefore, a 
state that maximizes the Hilbert-Schmidt norm of Y must be a pure state, 
Wr tS = 4>r4>* s , with \<ft r \ 2 = 1. Since \W r>s \ is symmetric with respect to 
exchanging r and s, we can symmetrize the term (1 — a r ) 2 a 2 in the above 
sum. 

To determine a maximizer, we use a variational principle. We replace the 
discrete spectrum of A with a continuous variable and examine the sym- 
metrized quantity 

1 1 
g(x, y) = xy(l - x)(l - y) + -(1 - x) 2 y 2 + -(1 - yfx 2 

= ^((1 - x)y + (1 - y)x) 2 . 
A critical point of g satisfies 

(l-2y)(y(l-x)+x(l-y))=0 

and the same equation with x and y exchanged. Thus, either x = y = or 
x = y = \. 

We denote k = \\A\\. 

Case I: k < 1/2. Considering the boundaries, we have g(x, 0) < g(n, 0) = 
|k 2 and re) = rr(l — x)k(1 — k) + |(1 — x) 2 re 2 + |(1 — k) 2 x 2 for < x < re, 
which is convex if 1 — 2re ^ and linear otherwise. In both cases the 
maximum on the boundary is either g(0, re) = |re 2 or g(re, re) = 2(1 — re) 2 re 2 . 
When re < |, the latter expression is the maximum. Moreover, the critical 
point is not contained in the open square (0, re) 2 . 
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This means that the state which maximizes 

tr[yy*] = 2 ^((1 - a r )a r (l - a s )a s + (1 - a r ) 2 a 2 )|0 r | 2 |0 s | 2 

r,s 

is given by being an eigenstate of A with highest eigenvalue k, and 

maxtr[yy*l = 4(1 - «)V . 
w 

Case II: 1/2 < k < 1. Since the rank of A is not maximal, we know it has an 
eigenvalue zero. From the preceding lemmas, we know that the probability 
vector p = 1 0| 2 which maximizes 

tr[yy*] = ^((1 - a r )a s + (1 - a s )a r ) 2 p r p s 

r,s 

can have at most two nonzero entries, which correspond to the maximal 
eigenvalue n and an eigenvalue zero. 

Choosing p accordingly, and optimizing 

tr[yy*] = (4p 2 (l - k) 2 k 2 + 2p(l - V )k 2 ) 

gives 

1 

P ~ 2-4(1 -k) 2 ' 

Inserting this, we obtain 

k 2 

maxtr[FF*] 



w 1 J 2-4(1-k) 2 ' 

After inspecting both cases for k, we conclude that 

n f A(\-k) 2 k 2 , ne [0,1/2] 
m - aXtr[FF] = ( «e (1/2,1] 



□ 



Theorem 3.4. Consider the Hilbert space C m <g)C' and /e£ £/ie set o/ general- 
ized phase-damping channels be the convex hull of the channels {£ n }™ =1 
defined with the set of projections {Q n = E nn ® I}™=i- F° r an U gi ven C*- 
encodmg $ : B(C d ) -> i3(C m <g> C'), Z < d ; tne worst case error zs bounded 
below by 

(16) sup e(*,£) > ( 2( ^^ )? 2 l [ ?:}' 2 }{ 

and equality holds if and only if all {V*Vj} are uniformly weighted projec- 
tions; that is, for all j G {1,2,... m}, we have V*Vj = kPj, with Pj = PjPj 
and k = 

ml 
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Proof. Given a C*-encoding $, we optimize among all convex combinations 
of S n as denned in the preceding lemma. 

We observe that for any probability vector p e K+, by Minkowski's in- 
equality for the Hilbert-Schmidt norm 

j j 

The right-hand side is maximized by the probability vector which is nonzero 
only in the index j with the largest e(<&,£j). However, choosing the worst- 
case state for this channel and the maximizing probability vector on the 
left-hand side gives equality in the above inequality. 

Therefore, optimal encoding suppresses the maximum among all e(<&, £„). 

The function e(<&, £„) is increasing in k n = \\A n \\. Since we want to choose 
{Aj} such that the maximum of e(<E»,£ n ) is minimized among all choices of 
n, we want to minimize maxi<j< m ||A,-||. The minimum is achieved when 
equality holds in Inequality (7). This requires that the POVM {Aj} consists 
of uniformly weighted projections. □ 

Remark 3.5. Recall that in Example 2.10 we saw how to construct uniformly 
weighted rank-one projections resolving the identity. Then in Example 2.11 
we indicated how to generalize this construction to higher-rank projections 
in the case that I divides d. We are aware of some other cases for which such 
families can be constructed, but here we leave the general case as an open 
problem that warrants further investigation. 



4. Optimality for next to minimal decoherence 

Next, we investigate what happens if we choose <£> optimal for minimal 
decoherence, and consider convex combinations of measurements with next 
to minimal impact. By the results of the preceding section, we can assume 
that the coordinate operators {Vj}™ =1 of $ factor a set of uniformly weighted 
projection operators, V*Vj = kPj, with k = d/ml and fcX^li Pj = I- 

Definition 4.1. Let K be a subset of two elements from the index set J = 
{1, 2, . . . m}. Define Q K = Y, je K E n ® J > and let f or M e B{®f =l Q), 

E K :M^ Q K MQ K + Q^MQj, . 

Let be the convex hull of all transmission errors £ K indexed by two- 
element subsets K of J . 

The following lemma reduces to the so-called Welch bound in the special 
case that all projections {Pj} are of rank 1 = 1. 
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Lemma 4.2. Let {P;}™ 1 be a set of rank-l projections forming a uniformly 
weighted resolution of the identity ^ Y^j=i Pj = I on a Hilbert space of 
dimension d > I. Then (a) if d < 21, we have that for all pairs i,j G 
{1,2, . . . m}, \\Pi + Pj\\ = 2; (b) if d> 21, on the other hand, then 



Im — d 



max WPi + Pj] > 1 + , . 
i±i V d[m - 1; 

and equality holds in this inequality if and only if for all i,j G {1,2,... m}, 



\Pi + Pj\\ = 1 + 



Im — d 
d(m — 1) 



Proof. We can assume that the Hilbert space is C d , equipped with the canon- 
ical basis, and identify the projections with d x d matrices. 

(a) We first show that in case d < 21, the sum of two of the projections has 
norm ||Pj + Pj\\ = 2 for all i,j G J. To see this, we write the projections as 
2x2 block matrices with the (1, l)-block being an I x I matrix, and choose 

a basis such that this block is the identity for Pi = ■ Now when 

we compress P 3 to this first I x / block, we remove r = d — I < I rows and 
columns and so by eigenvalue interlacing, the / x / block still has eigenvalue 
1 of multiplicity I — r > 0. Consequently, Pj + P, has 2 as an eigenvalue. 

(b) We now consider the case d > 21. 

Splitting C d = C l © © C r , r > 0, we write dx d matrices as 3 x 3 block 
matrices. By a suitable choice of basis, we may assume that 

// N 
Pi = 
\0 0, 

In the same decomposition, each Pj is written as a block matrix, and we call 
Bj the (1, l)-block entry, so Bj is an I x / matrix. 

Now summing these matrices we have that / + Y^j=2 Bj = if I by ^ e 
resolution of the identity and hence X]j=2 tr (-^j) = l ~f~ ~ ^ Thus, selecting 
for each B 3 its largest eigenvalue \ 3 gives the chain of inequalities 

Im — d 



(17) maxA, > A, > 

j>2 J ~ m - 1 4^ ~ 



d(m - 1) • 

In the remaining part of the proof, we show that 

||Pl + P,|| = 1 + y/Xj. 
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Once this is established, we have by the monotonicity of the square root that 



Im — d 



max 1 1 Pi + Pj\\ > 1 + , ( , 
i>2 J y d(m - 1) 

and equality implies that equality holds in the chain of Inequalities (17), 
which in turn gives that all Xj are equal, as well as all ||Pi + P,||. Since the 
choice of Pi is arbitrary, this argument shows that all pairs {P + Pj}i^j have 
equal norm. 

For simplicity of notation, we continue the remaining part of the proof 
with i = 1 and j = 2. We first show that up to unitary equivalence P 2 = 
' D R \ 

R R 2 2 R23 , with diagonal blocks D, R, and as yet unspecified parts 
P 32 P33 J 

R22, R23, R32, and P33. To see this, first we write P 2 in 2 x 2 block form with 
the (1, l)-block an / x / matrix and conjugate by a unitary to diagonalize B 2 
as D. 

Now we polar decompose the (1, 2)-block as RV* where R > is / x I and 
V* is I x (d — I), and by assumption d — I > I. Let U be a (d — I) x (d — I) 

unitary so that V*U = (10). Conjugating by ^ j T ) we have that up to 

/ D R 

these two unitary transformations, P 2 = I R R22 R23 \ ■ But since P 2 is a 

\0 P32 

projection, D 2 + R 2 = D, and we see that R is diagonal as well. 

If R has entries that are zero, then the corresponding entries in D must 
vanish also, since by assumption ||P! + P 2 || < 2 which yields \\D\\ < 1. 
Again using that P 2 is a projection, DR + PP 22 = R and hence P 22 is also 
diagonal in the rows and columns where R is non-zero. In addition, since 
= RR23, R23 vanishes in those rows. 

We now observe, similarly as in Example 2.12, that P 2 is a direct sum of 
2x2 matrices and a piece which contains the blocks that are not diagonal. 

If the eigenvector of Pi + P 2 corresponding to the largest eigenvalue had 
nonzero entries in a row in which D and R vanish, then that eigenvalue 
would be one. Indeed, this is the case if P 2 = D = R = 0, since then the 
ranges of Pi and P 2 are orthogonal, and ||Pi + P 2 || = 1. 

However, if D has rows with non-zero entries, then we find the largest 
eigenvalue among the 2x2 matrices. To this end, pick a non-zero entry 
of D, say A. Note that since P 2 is a projection, the corresponding diagonal 
entry of R is then y/X — X 2 and that of P 22 is 1 — A. Thus, in P 1 + P 2 we have 

a 2 x 2 matrix of the form f ~\2 ^ ^ ' ^ ne finds that the norm 
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of this matrix is 1 + y/\. Consequently, the largest entry A2 in D maximizes 
this norm, which shows that ||Pi + P 2 \\ = 1 + || V^B^H- D 

Theorem 4.3. Assume the C* -encoding $ : B(C d ) -> B(C m <g> C') has 
coordinate operators {V 7 }JL 1 with rank I < d, such that each V*Vj = kPj 
with a projection operator Pj and k = d/ml. Denote 



k = 



d 1 / d(lm-d)~ 1 - 9? 
ml + V m^(m-l) ' " - U ' 



£/ien £/ie worst case error under the C* -encoding $ and transmission channels 
from Q( 2 ) is bounded below by 

,~ ^ f 2«(l-«), KG [0,1/2] 



If d > 21, then the value max £ - 6 g( 2 ) e($, 5) attains the lower bound if and 
only if for all i,j G {1,2,... m}, 



\Pi + Pi\\ = 1 + 



Im — d 



dim 



Proof. Denote A = kYlj^K^j^i- After decoding, the resulting error in the 
reconstructed state is 

Y = (I - A)WA + AW(I - A) . 

In the eigenbasis of the operator A, the square of the Hilbert-Schmidt norm 
of Vis 

tr[yy*] = 2 ((! - «r)«r(l - <x a )a 8 \W r>9 \ 2 + (1 - a r ) 2 «^r,.| 2 ) , 

where again the matrix with entries (W rjS ) represents the state W in the 
eigenbasis of A and {o; r } are the corresponding eigenvalues of A. 

We notice that finding the state W which maximizes the worst-case error 
norm is exactly the same problem as in Lemma 3.3 of the preceding section. 
By the monotonicity of the resulting expression in \\A\\ and the convexity 
of the Hilbert-Schmidt norm, we can conclude analogously that the optimal 
encoding minimizes the largest of the eigenvalues of all pairs {Pi + Pj}i^j. 
The preceding lemma states that this is the case if and only if the norms of 
all these pairs are equal. The common value of their norms then yields the 
claimed error bound. □ 



We describe a class of examples for which the lower bound for the worst- 
case error is attained. 
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Example 4.4. Let C d = C q © C l , and {Pj = Uj © I}f =1 with rank-one 
projections {LL,} giving a uniformly weighted resolution of the identity on 
the first component C g . We follow ideas of Holmes and Paulsen [25] to 
characterize optimal uniformly weighted projections of this type. 

Let {fjjJLi be a set of non-zero vectors, such that each fj is contained in 
the span of the corresponding rank-one projection H,, with the normalization 
WfjW = V q/m. One can check that this set is a Parseval frame. 

The largest eigenvalue of Pi+Pj is + = ||(rij+rij)(g)/|| = ||H; + n,-||. 
From identities for the trace of the Grammian G of the Parseval frame {fj}, 
we have q = tiG = E™i \\fj\\ 2 and q = trG*G = EijKfjJ^ 2 , and we 
can infer the usual form of the Welch bound [25,27,36] 

max ||n, + n.H = max {1 + -\(f, fj)\} > 1 + 
i+3 q 

and equality holds if and only if \(fi, fj)\ = \J q(m — q)/ (m 2 (m — 1)) for all 

Therefore, an optimal set of projections of the form {LL; <g> 1} with rank- 
one projections {LL,} is given by a so-called two-uniform frame {fj} with all 
inner products between two different frame vectors having the same absolute 
value. 

Remark 4.5. Now we briefly remark how to generalize to the case when 
dimV = s > 0. This includes the case m — 1 of so-called decoherence- 
free subspaces. The minimal error model is the convex hull of transmission 
channels given by projections of both types {Qj = (Ejj ® I) © 0}JLi or 

{Qm+n = (0 <g> 0) © E nn } s n=l . 

Suppose we have a C*-encoding $ : B(C d ) -> B((C m © C') © V). Then we 
can view identify the tensor product as a direct sum, (C m © C l ) © V = 
{(BjLiC 1 ) © V and split the matrix V associated with $ similarly as in 
the proof of Lemma 2.6 into coordinate operators {Vj}™^ of rank / and 
{Vj}fJ^ +l of rank one. 

We note that the coordinate operators give rise to a POVM {Aj}™^ 1 with 
m rank-/ operators {Aj}™ =1 and s rank-one operators {^4j}j=tn + i- 

In analogy with Theorem 3.4, we can then derive a modified version of 
the lower bound for the worst case error, which amounts to replacing ml 
with ml + s in Ineq. (16). Again, equality would hold for $ with coordinate 
operators that factor uniformly weighted projections. However, apart from 
the case when s divides I it is not clear under which general conditions an 
optimal $ can be constructed. 

Similarly, to find the best $ for the next to minimal error model, we 
require the maximal eigenvalue among all {Pi + Pj}i^j to be minimized, but 



m — q 
q(m — 1) 
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at this point it is not clear whether there are enough interesting examples of 
this type. 
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